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Abstract: The Inverted Pendulum is a very popular plant for testing dynamics and control of highly non-linear 
plants. In the Inverted Pendulum Control problem, the aim is to move the cart to the desired position and to 
balance a pendulum at desired location. This paper represents stabilization of pendulum using PID, SVFB and 
IQR. An advantage of Quadratic Control method over the pole-placement techniques is that the former 
provides a systematic way of computing the state feedback control gain matrix.IQR In the mathematical model 
proposed here, a PID, SVFB and IQR controller is designed and all the controllers are tuned using Genetic 
algorithms. The simulation results of all the controllers are shown. This paper exhibits the capability of genetic 
algorithm to solve complex and constrained optimization problems and as a general purpose optimization tool 
to solve control system design problems. 

Keyword: Inverted pendulum (IP), Mathematical modeling, PID controller, State variable feedback controller 
(SVFB), IQR, Genetic algorithm (GA). 



I. INTRODUCTION 

The inverted pendulum is unstable [8,9,10] in the sense that it may fall any time in any direction unless a 
suitable control force is applied. If the designer works it right, he can get the advantages of several effects [11]. 
The control objective of the inverted pendulum is to swing up[4] the pendulum hinged on the moving cart by a 
linear motor[ll] from stable position (vertically down state) to the zero state (vertically upward state) and to 
keep the pendulum in vertically upward state in spite of the disturbance[6,7]. They make it easy to check whether 
a particular algorithm [5] yields the requisite results. Several works has been reported on the inverted pendulum 
for its stabilization. Attempts have been made in the past to control it using classical control [3]. The purpose of 
the present research is to tune[l,2] the entire controller by soft computing tools. The work was made under 
MATLAB simulation. 



II. MATHEMATICAL MODEL OF THE PLANT 

Defining displacement of the cart as x, the angle of the rod from the vertical (reference) line as 9 , the 
force applied to the system as F , centre of gravity of the pendulum rod at its geometric centre and I the half 
length of the pendulum rod, the physical model of the system is shown in fig (1). 
The Lagrangian of the entire system is given as, 

L= \prtx 1 +2mlx9cos9+ml 2 9 2 +Mx 2 )+ ^I9 2 ]-mglcos9 

The Euler-Lagrange's equation for the system is 

±('±)-'± + bi = F (i) 

dt \SxJ Sx y ' 

±( s ±)_ s ± +d g = 0 (2 ) 
dt \seJ se v ' 

The dynamics of the entire system using above equation is 

(/ + ml 2 )9 + mlcos9x- mgl sin 9 + d9 = 0 (3) 
(M +m)x + ml cos 9 9 -mlsm9 9 2 + bx = F (4) 

In order to derive the linear differential equation model, the non linear differential equation obtained need to be 
linearized. For small angle deviation around the upright equilibrium (fig. 2) point, assumption made sin 9 = 
9, cos 9 = 1, 9 2 = 0 

Using above relation, equation (5) and (6) are derived. 

r0'+qx-k0+d§=O (5) 

px +q9 + bx = F (6) Where, (M 

+ m)= p, mgl=k, ml=q, I + ml 2 = r 
x 1 — x + / sin 9 
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Fig 1 : The Inverted Pendulum System 
Table 1 . Parameters of the system from feedback instrument .U.K. 



Parameter 


Value 


unit 


Cart mass(M) 


1.206 


Kilo gram 


Mass of the 


0.2693 


Kilo gram 


pendulum(m) 






Half Length of 


0.1623 


meter 


pendulum(Z) 






Coefficient of 


0.005 


Ns/m 


frictional force(fo) 






Pendulum 


0.005 


Mm/rad 


damping 
coefficient(d) 






Moment of inertia 


0.099 


kg/m 2 


of pendulum(Z) 






Gravitation 


9.8 


mis 2 


forced) 







Eq (5 & 6) is the linear differential equation model of the entire system. 

After taking Laplace transform and substituting the parameter value from (table 1), we got 

8(s) -0.2783 S 2 



F(s) 
X(s) 

F(s) s(s+2.026)(s 



s (s+2.026) (s-1.978) (s+0.03402) 
0.68843 (s+2.014) (s-1.967) 



: G,(S) 

: G 2 (S) 



1.978) (s+0.03402) 

3 STATE SPACE MODELLING 

Lct ? 0=A^ ? Q—x-^—x-^^ Q—x-^—'x-^ iuid, x — x-$, x—x^—x-^^ x—x^—x-^ 

From state space modeling, the system matrices are found in matrix form as given below. 
0 1 0 0.0000000' 

4.0088 -0.047753952 0 0.0139155 



(7) 
(8) 



0 



0 



0.116806166 0.0939155 



0 
0 



1.0000000 
0.0344200 



and 



B= 



0 

-0.27831 
0 

0.68842 



rl 0 0 0] 

"1.0 o i oJ 



*3 



Y be the output equation 



III. PID Controller design 

PID controllers are a family of controllers. The reason PID controllers are so popular is that using PID gives the 
designer a large number of options and those options mean that there are more possibilities for changing the 
dynamics of the system in a way that helps the designer. 
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*-r->- 



n= Output forcing function = Kp-C +Kj-^,+ K||*dt 



H[S)=1 



Fig 2: Design of PID controller 
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Fig 3: Design of PID controller 



IV. LQR DESIGN 

A system can be expressed in state variable form as 
x = Ax + Bu 

The initial condition is x(0) . We assume here that all the states are measurable and seek to find a state-variable 
feedback (SVFB) control 
u = —Kx + v 

that gives desired closed-loop properties. 

The closed-loop system using this control becomes 

x = (A- BK) x+ Bv= A c x+ Bv 

With A c the closed-loop plant matrix and v the new command input. 

Ackermann's formula gives a SVFB K that places the poles of the closed-loop system at desired locations. To 
design a SVFB that is optimal, we may define the performance index J as 

J = 1 -fiX T (Q+K T RK)Xdt 

We assume that input v(t) is equal to zero since our only concern here are the internal stability properties of the 
closed loop system.The objective in optimal design is to select the SVFB K that minimizes the performance 
index J. The two matrices Q and R are of appropriate dimension. One should select Q to be positive semi- 
definite and R to be positive definite. Since the plant is linear and the performance index is quadratic, the 
problem of determining the SVFB K to minimize J is called the Linear Quadratic Regulator (LQR). To find the 
optimal feedback gain matrix K we proceed as follows. Suppose there exists a constant matrix P such that 

— (x T Px) = —x T (Q + K T RK)X 
at 



After some mathematical manipulation, the equation becomes, 
1 C x d 1 

y = - 2 J ^(* r ^) dt = 2* r(0)p * (0) 



Where, we assumed that the closed-loop system is stable so that X(t) goes to zero as time t goes to infinity. 
Substituting the values we get, 
x T (A T c P + PA C + Q+K T RKx)x = 0 

It has been assumed that the external control v(t) is equal to zero. Now note that the last equation has to hold for 
every X T . Therefore, the term in brackets must be identically equal to zero. Thus, proceeding one sees that 



www.ijeijournal.com 



Page I 70 



Tuning of PID, SVFB and LQ Controllers Using Genetic Algorithms 



A T P + PA + Q+ K T RK - K T B T P - PBK = 0 

Exactly as for the scalar case, one may complete the squares. Though this procedure is a bit complicated for 
matrices, suppose we select 
K = R~ 1 B T P 
Then, this results with 
A T P + PA + Q- PBR~B T P = 0 
This is algebraic matrix Riccati equation 

V. GENETIC ALGORITHM 

Genetic algorithms (GA) are search procedures inspired by the laws of natural selection and genetics. They can 
be viewed as a general-purpose optimization method and have been successfully applied to search, optimization 
and machine learning tasks. GA has the ability to solve difficult, multi dimensional problems with little 
problem-specific information and hence has been chosen as the optimization technique to solve various 
problems in control systems.lt has been shown that compared with other traditional heuristic optimization 
method, Genetic Algorithm is likely to be more computationally efficient. The controller parameters are usually 
determined by trial-and-error through simulation. In such case, the paradigms of GA appear to offer an effective 
way for automatically and efficiently searching for a set of control performance. 



The flow chart for implementing GA for tuning controller parameters is shown below. 



Create/Initialize 
Population 



Measar&'Evaluate 

Fitness 



Crossover ■' Production 



Non optimum 
Solution 



Optimum Solution 



Fig 4:Flow chart of Genetic Algorithm 

The steps involved in creating and implementing a genetic algorithm: 

(a) Generate an initial, random population of individuals for a fixed size. 

(b) Evaluate their fitness. 

(c) Select the fittest members of the population. 

(d) Reproduce using a probabilistic method (e.g., roulette wheel). 

(e) Implement crossover operation on the reproduced chromosomes 

(f) Execute mutation operation with low probability. 

(g) Repeat step 2 until a predefined convergence criterion is met. 



VI. SIMULATION AND RESULTS 



6.1 GA based PID controller: 




Fig 5 GA based tuning with and without of PID Controller 
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Response of plant with PID Controller 
Response of plant with GA tuning 



Results:- For 100 Population 

Kp = 50.8689 Kd = -58.9523 Ki = -149.2234 

6.2 Genetic Algorithm based State feedback Controller: 




Response of IP System 

Response of State Feedback 



Controller with GA tuning 



0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



Fig 6 GA based tuning with and without of State Feedback Controller 



6.3 Genetic Algorithm based LQR Controller 

-3 

X 10 




Response of IP System 

Response of IP with GA tuning 



10 12 14 16 18 20 



Figure 7 GA based tuning with and without of LQR 

VII. CONCLUSION 

Modeling of inverted pendulum shows that system is unstable with non-minimum phase zero. Results 
of applying PID controllers show that the system can be stabilized. While PID controller method is cumbersome 
because of selection of constants of controller, Constant of the controllers can be tuned by some Genetic 
Algorithm technique for better result. Results with GA tuned controller are better in respect of rise time and 
overshoot.Fuzzy Logic Controllre may be used for best result. 
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